Astronomy 


& Astrophysics manuscript no. MS printer 


© ESO 2009 


November S 


, 2009 





Photochemistry of the PAH pyrene in water ice: the case for 
ion-mediated solid-state astrochemistry 

' Jordy Bouwman 1 , Herma M. Cuppen 1 , Arthur Bakker 1 , Louis J. Allamandola 2 , and Harold Linnartz 1 

o ■ 

0\ 



Raymond and Beverly Sackler Laboratory for Astrophysics, Leiden Observatory, Leiden University, P.O. Box 9513, NL 2300 RA Leiden, 
The Netherlands 

2 NASA-Ames Research Center, Space Science Division, Mail Stop 245-6, Moffett Field, CA 94035 
Received ; accepted 

ABSTRACT 



on 

: 

I • Context. Icy dust grains play an important role in the formation of complex inter- and circumstellar molecules. Over the past decades, laboratory 
O ' studies have mainly focussed on the physical interactions and chemical pathways in ices containing rather simple molecules, such as H 2 0, CO, 
CO2, CH 4 and CH 3 OH. Observational studies show that polycyclic aromatic hydrocarbons (PAHs) are also abundantly present in the ISM in the 
gas phase. It is likely that these non-volatile species freeze out onto dust grains as well and participate in the astrochemical solid-state network, 
but experimental PAH ice studies are largely lacking. 

Aims. The study presented here focuses on a rather small PAH, pyrene (Ci 6 H 10 ), and aims to understand and quantify photochemical reactions 
of PAHs in interstellar ices upon vacuum ultraviolet (VUV) irradiation as a function of astronomically relevant parameters. 
Methods. Near UV/VIS spectroscopy is used to track the in situ VUV driven photochemistry of pyrene containing ices at temperatures ranging 
from 10 to 125 K. 

Results. The main photoproducts of VUV photolyzed pyrene ices are spectroscopically identified and their band positions are listed for two 
host ices, H 2 and CO. Pyrene ionisation is found to be most efficient in H 2 ices at low temperatures. The reaction products, triplet pyrene 
and the 1 -hydro- 1-pyrenyl radical are most efficiently formed in higher temperature water ices and in low temperature CO ice. Formation routes 
and band strength information of the identified species are discussed. Additionally, the oscillator strengths of Py, Py + and PyH are derived and 
a quantitative kinetic analysis is performed by fitting a chemical reaction network to the experimental data. 

Conclusions. Pyrene is efficiently ionised in water ice at temperatures below 50 K. Hydrogenation reactions dominate the chemistry in low 
temperature CO ice with trace amounts of water. The results are put in an astrophysical context by determining the importance of PAH ionisation 
. in a molecular cloud. We conclude that the rate of pyrene ionisation in water ice mantles is comparable to the rate of photodesorption of H 2 
' ice. The photoprocessing of a sample PAH in ice described in this manuscript indicates that PAH photoprocessing in the solid state should also 
, be taken into account in astrochemical models. 

Key words. Spectroscopy: Solid State, methods: laboratory 



C3 



> 

O 
in 



On 
O 



1. Introduction 

Strong infrared emission attributed to polycyclic aromatic 
hydrocarbons (PAHs) is character i stic o f many galactic and 



extragalactic objec ts (Smith et al. 2007t iDra 
Draine & Lil |2007t iTielensI I200C While this 



Draine et al. 2007 



emission gen- 
erally originates from optically thin, diffuse regions, PAHs 
should also be common throughout the dense interstellar 
medium. There, as with most other interstellar species in 
molecular clouds, PAHs condense out of the gas onto cold icy 
grain mantles where they are expected to influence or partici- 
pate in the chemistry and physics of the ice. While laboratory 
studies on interstellar ice analogs have shown that complex or- 
ganic molecules are produced upon extended vacuum ultravi- 
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olet ( VUV) photolysis (e.g. lBriggs et al.lll992tlBernstein et al 
1995b . the photoinduced processes at play during the irradia- 
tion of PAH containing interstellar ice analogues have not been 
studied in detail. In optical, in situ studies on the photochem- 
istry of naphthalene, 4-methylpyrene and quatterylene contai n- 
ing water ice at 20 K , iGudipati & Allamandolal d2003l l2006albl) 
and iGudipatil (120041) showed that these PAHs are readily ion- 
ized and stabilized within the ice, suggesting that trapped ions 
may play important, but overlooked roles in cosmic ice pro- 
cesses. Beyond this, there is little information on the VUV in- 
duced, in situ photochemistry and photophysics of PAH con- 
taining water rich ices. 

Here, we describe a detailed study of the VUV in- 
duced photochemistry that takes place within pyrene (Py or 
Ci 6 Hio) containing water ices (Py:H 2 O=l:10,000-l:5,000). 
The present study is an extension of a recently published study 
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(Bou wman et al. 2009) in which the focus was on the new ex- 
perimental setup and where the use of PAH ice spectra was 
discussed to search for solid-state features of PAHs in space. 
In this work, the focus is on a detailed characterization of 
the actual chemical processes taking place upon VUV irra- 
diation, particularly as a function of ice temperature ranging 
from 25 to 125 K. Additionally, measurements on Py:CO ices 
at 10 K have been performed to further elucidate the role of 
water in the reaction schemes and to clarify formation routes 
of identified species. A similar study of three small PAHs 
is now underway to understand the general principles of 
PAH/ice photochemistry. This is part of an overall experi- 
mental program at the Sackler Laboratory for Astrophysics 
to study the fundamental processe s of inter- and circum stel- 
lar ice analogues such as thermal dAchary va et al.l 120071) and 
photodesorption ( Ober g et al1l2007bl l2009h. hydogenat ion re- 
actions (IFuchs et alJIacce pted 2009t lloppolo et al.ll2008l) . pho- 
tochemistry dOberg et al.l 12.009b and physical in t eractions in 
interstellar ice analogues (Bouwman et al . 2007; Ober g et al 
2007ail2009l) . 



phot on flux of ~ 10 14 photons-cm 2 s at the sample surface 



The manuscript is organized as follows. The experimen- 
tal technique is summarized in Section [2] Section [3] describes 
the Py:H20 and Py:CO ice photochemistry, the resulting prod- 
ucts and their formation routes. The temperature dependent 
photochemistry and derived reaction dynamics are described 
in Section [4] and astrochemical implications are discussed in 
Section[3] The main conclusions are summarized in Section[6] 



2. Experimental Technique 

We use a recently constructed apparatus, described in 



Bouwman et al.l ( 20091) . that follows the photochemistry in ki- 
netic mode during VUV irradiation by measuring the near 
UV- Visible absorption spectra of an ice, providing 'real-time' 
tracking of the reactants and photoproducts. Dilute Py:H20 ice 
samples (-1:10,000 - -1:5,000) and a Py:CO ice sample of 
comparable concentration are prepared by depositing the vapor 
from a pyrene sample heated to 40°C together with H2O vapor 
or CO gas onto a cold MgF2 window. The window is cooled to 
10 K in the case of CO deposition or 25 K in the case of H2O 
deposition. The sample window is cooled by a closed cycle He 
refrigerator. Pyrene (Aldrich, 99%) and CO (Praxair 99.999%) 
are used without further purification. Vapor from water, filtered 
through a milli-Q purification system and further purified by 
three freeze-pump-thaw cycles, is used. The sample window 
is mounted in a high-vacuum chamber (P a 10~ 7 mbar). Ice 
growth rate and thickness are determined with a HeNe laser 
by monitoring the thin-film interference fringes generated dur- 
ing deposition. Simultaneously, the amount of pyrene is tracked 
measuring the integrated strength of the S2 <— So neutral Py 
transition at 334 nm. Deposition is typically stopped when the 
optical depth (OD) of Py approaches ~0. 15. 

The ice samples are photolyzed with the 121.6 nm Lya 
(10.6 eV) and the 160 nm molecular hydrogen emission bands 
(centered around 7.8 eV) generated by a microwave powered 
discharge in a flowi ng H2 gas with a vacuum ultraviolet flux of 
~ 10 15 photons s _1 dMunoz Caro et alj|2002l) . This results in a 



dOberg etaDl2009h . 

Absorption spectra of VUV photolyzed Py containing ices 
are measured with a Xe-arc lamp serving as a white light 
source. Lenses and diaphragms direct the light through the ice 
sample along the optical axis determined by the HeNe laser 
beam after which it is focused onto the entrance slit of a 0.3 m 
spectrometer. A 150 lines mirT 1 grating, blazed at 300 nm, dis- 
perses the light onto a sensitive 1024x256 pixel CCD camera 
with 16 bit digitization. The camera is read out in the vertical 
binning mode by a data acquisition computer on which the data 
are converted to absorbance spectra (OD = -ln(///o))- This 
configuration spans the 270 to 830 nm spectral range, which 
permits simultaneously monitoring the behavior of the neutral 
Py parent molecule and photoproduct bands without adjusting 
any element along the optical path. This is critical to obtain reli- 
able and reproducible baselines in measuring the optical spec- 
tra of ices. The spectral resolution is of the order of 0.9 nm 
which is more than sufficient to record broad solid-state ab- 
sorption features. 

The measurements described here were made on various 
H 2 0:Py ice samples at 25, 50, 75, 100 and 125 K. The CO ice 
experiments were carried out at 10 K to avoid matrix sublima- 
tion at higher temperatures. The sample temperature is main- 
tained using a resistive heater with an accuracy of +2 K. The 
measured spectra are converted into units of optical depth us- 
ing the spectrum of the freshly deposited, unphotolyzed ice 
at the appropriate temperature as a reference spectrum (/o). 
Recording a single spectrum typically takes about 5 ms and 229 
spectra are generally co-added to improve the S/N of a spec- 
trum, resulting in one complete spectrum every 10 seconds. 

The optical configuration of the apparatus is such that 
spectra are recorded simultaneously with photolysis. Thus, the 
fast spectral recording time permits monitoring photoinduced 
changes on a roughly 10 second time scale. Figure Q] shows the 
290 to 490 nm spectrum of a Py:H 2 ice at 125 K after 900 s 
of in situ VUV photolysis. Because the spectrum recorded be- 
fore VUV irradiation is taken as a reference (Iq), bands with 
positive OD values originate from species that are produced by 
photolysis while the bands with negative OD correspond to the 
neutral pyrene that is lost upon photolysis. Comparing the Py 
and photoproduct absorption bands with the narrow H/3 lamp 
line at 486.1 nm shows that the instrumental resolution indeed 
far exceeds the ice band widths. The absolute wavelength cali- 
bration is accurate to within ±0.5 nm. 

More than 1400 individual spectra are recorded and are re- 
duced in a typical 4 hr experiment. Spectra are individually 
baseline corrected by fitting a second order polynomial through 
data points where no absorptions occur and subsequently sub- 
tracting the fit from the measured spectrum. Integrated ab- 
sorbances of absorption features are calculated numerically for 
all spectra. These are corrected for contributions by atomic 
hydrogen lines originating from the H2 discharge lamp. The 
data reduction software also allows plotting of correlation di- 
agrams between integrated absorbances of different absorption 
features. All data handling and reduction is done with Lab View 
routines. 
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Fig. 1. The spectrum of a dilute pyrene:H20 ice after 900 s of VUV irradiation at 125 K. The inset shows a blow-up of the pyrene 
photoproduct bands. Band assignments are discussed in Section [3] Note the broad feature ranging from about 350 to 470 nm 
which is indicated by a Gaussian fit. This is attributed to overlapping bands from individual pyrene photoproducts. Bands with 
negative optical depth indicate species destruction, those with positive optical depth show species formation. The blue bands are 
Gaussian profiles which co-add to the overall fit shown in red. Note the instrumental resolution indicated by the profile of the H/? 
line at 486. 1 nm. (This figure is available in color in electronic form) 



Integrated band areas are used, in conjunction with os- 
cillator strengths (/), to derive molecular abundances. The 
oscillator strength is converted to integrated absorbance 
(cm molecule" 1 ) usin g the conversion factor 8.88xlO -13 
(IKjaergaard et al . 2000). The number of molecules per cm 2 (AO 
is given by: 



N = 



f V2 r dv 



8.88 x 10- 13 /' 



(1) 



where t is the optical depth and v is the frequency in cm 



3. Band assignments and band strength analysis 

The typical photolysis duration of about 4 hours is the time re- 
quired for nearly complete loss of the neutral pyrene vibrational 
progression at 334.0, 329.2 and 319.2 nm. Irradiating the sam- 
ple ices with VUV light produces a set of new absorption bands 
in the spectra, indicating active photochemistry. The band posi- 
tions, FWHM and assignments of the bands in the Py:H20 ice 
at 25 K are listed in TableQ] The bands appearing in the Py:CO 
ice at 10 K are similar to those in the Py:H20 ice at 25 K, how- 
ever, with slightly altered band positions and FWHM and with 
very different relative intensities (see also Table [TJ. Figure Q] 
presents a spectrum from the 125 K Py:H20 series. This figure 
illustrates production of the pyrene radical cation (Py + ), triplet 
pyrene ( 3 Py), 1 -hydro- 1 -pyre nyl radical (PyH) and a broad un- 
derlying 'residue' feature upon VUV irradiation. Additionally, 
a progression of distinct absorptions is found in the Py:CO ex- 
periment, which indicates the formation of the (reactive inter- 
mediate) HCO radical. The identifications of these species and 
their oscillator strengths are discussed below. 



Table 1. Band positions (A c ) and FWHM of Gaussian fits to the 
observed features of neutral pyrene for pure pyrene ice at 10 K, 
pyrene in a water ice at 25 K and pyrene in a CO ice at 10 K. 
Photoproduct bands arising from VUV processing of the ices 
are also listed for the Py:H20 and Py:CO experiments. Both 
band positions and FWHMs are indicated in nm. 



Species 


Pyrene 


Pyrene:H20 


Pyrene:CO 




A c FWHM 




FWHM 


A c 


FWHM 


Py 1 B 2 „ 


312.7 7.1 


319.2 


6.5 


319.4 


7.5 




325.3 10.0 


329.2 


3.2 


329.2 


2.3 




341.5 14.0 


334.0 


4.4 


334.3 


4.1 


Py' + 2 B hl 




363.2 


3.6 










354.0 


6.5 










344.9 


6.2 






Py + 2 A„ 




445.6 


6.6 


445.3 


7.8 






435.5 


10.2 










423.0 


12.2 










413.8 


5.3 






Py + 2 B h , 




490.1 


10.0 






PyH 




399.4* 


5.2'' 


400.5 


4.2 










392.5 


6.7 










378.4 


15.7 


3 Py 3 A- 




405.0* 


4.5* 


406.2 


4.8 


HCO i A" 








513.4 


17.5 










535.3 


12.5 










556.3 


14.5 










583.0 


16.8 










604.9 


10.0 










639.2 


15.1 



17 Absorption feature was too weak to perform an accurate fit 
Features are too weak at 25 K; the 125 K values are indicated 



3.1. Neutral pyrene bands 

As in Bouwman et al. (2009), the strong, negative bands peak- 
ing at 334.0 nm and weaker bands at 329.2 and 319.2 nm in 
the H2O ice (see Fig.[TJ, and at slightly shifted positions in the 
CO ice, are assigned to the 1 B2« <_ l A g electronic transition of 
neutral pyrene (S2 <— So) based on previous studies of pyrene 
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in rare gas matrices dVala et al.lll994t lHalasinski et alj|2005l) . 
To study the chemistry in absolute number densities, a value 
of / = 0.33 is ado pted from the li t erature for the osc illator 
strength of pyrene jBito et all [20001: IWang et al.ll2003l) . This 
value is used throughout this paper both for the Py:H20 and 
Py:CO experiments. Pure pyrene ice measured at 10 K shows 
broader absorptions located at 341.5, 325.3, and 312.7 nm (see 
Table [T]i. We did not perform VUV experiments on the pure 
pyrene sample. 

3.2. Pyrene cation bands 

Positive bands at 363.2, 354.0 and 344.9 nm appear upon pho- 
tolysis in the Py:H20 experiments. This progression is assigned 
to the 2 Bi„«— 2 Z?3 g vibronic transition of the pyrene cation (Py + ) 
in accordance with th e pro ximity to the band posi tions reported 
by IValaetalJ dl994l) and lHalasinski etaD J2005I). This tran 



sition for Py' + in H2O ice was reported in iBouwman et al 
d2009l) . The 2 .Bi u <— 2 B3 ? transition is too weak to be detected 



in the Py:CO experiment. A stronger Py + progression falls at 
445.6, 435.5, 423.0, and 413.8 nm in water ice. Of these bands 
only the strongest at 445.3 nm is detectable in the irradiated 
Py:CO ice. This progression is assigned to the 2 A„<— 2 B^ g of 
Py' + . The much weaker absorption due to the 2 Z?i„<— 2 Z?3 g Py' + 
transition at 490.1 nm in H2O is again undetectable in CO. 

Py + formation in these H2O and CO ice experiments is the 
result of direct single photon ionisation of the neutral species, 
following: 



Py 



vuv 



Py + + e" 



(2) 



It has to be emphasized that ionisation in Py:H20 ices is far 
more efficient than in Py:CO ices. Additional measurements 
on Py:CO:H20 mixtures indicate that the presence of H2O in- 
deed enhances the ionisation. Hence, it is possible that water 
contamination in the CO ice is responsible for the formation of 
some, if not all, of the cation species in the Py:CO experiment. 
The role of water contamination in CO ice will be discussed in 
more detail in Section [3~3l 

Using baseline corrected spectra as shown in Fig. [T] the 
photochemical evolution is tracked by integrating areas of 
bands due to each species in every spectrum and plotting them 
as a function of photolysis time. The strongest Py' + band at 
~445 nm is selected to track the number density evolution 
of this species. To put the kinetic analysis (Section HI on a 
quantitative footing, we determine the oscillator strength of 
the 445 nm Py + band as follows. First the integrated ab- 
sorbance of the 445 nm Py' + band is plotted versus that of the 
334 nm Py band during the course of VUV photolysis at differ- 
ent ice temperatures. These graphs are shown in Fig. [2] Note 
the tight, linear behavior between the loss of neutral pyrene 
and growth of the pyrene cation during early photolysis times 
up to 100 s (the first 10 successive datapoints). Inspection of 
Fig. [2] shows that the slope is steepest and the ratio of the in- 
tegrated absorbance of the Py' + band to the Py band is opti- 
mum in the 25 K ice. Since no other photoproduct bands are 
evident during the linear correlation stage, we conclude that 
during this phase, neutral pyrene is solely converted to the 
cation as described previously for naphthalene and quaterry- 




-0.1 -0.2 -0.3 -0.4 -0.5 -0.6 
Normalized integrated absorbance Py neutral 

Fig. 2. Integrated absorbance of the 445.6 nm Py' + band growth 
plotted against the loss of the 334.0 nm Py band in the 25, 50, 
75, 100, and 125 K ices. 10 seconds elapse between subsequent 
data points. Values are normalized to the maximum integrated 
absorbance of neutral pyrene. The straight line portion of these 
plots is used to determine the oscillator strength of Py + as de- 
scribed in Section l3~2l 



lene (Gudipat i & Allamandolal 2006a ). The straight line por- 
tion, fitted through the first 10 data points of irradiation at 25 K, 
is used to determine the oscillator strength of Py + . Given that 
the ratio of the Py + to the 334 nm Py band is 0.99 and the 
oscillator strength of this Py transition is 0.33, the oscillator 
strength of the 443 nm Py' + band in water ice is taken to be 
0.33 as well. This conclusion is consistent with ab initio cal- 



culations on pyrene bv lWeisman et al.l (120051) . They calculated 
that the oscillator strength of the cation is only ~2% stronger 
than that of the neutral species. 

As will be described below, the photolysis of Py in water 
ices at higher temperatures produces other species in addition 
to the cation. This explains the different curves in Fig. [2] 

3.3. HCO bands in Py.CO 

VUV irradiation of a Py:CO ice also produces a vibrational 
progression ranging from ~ 500 to 650 nm. As shown in Fig.[3j 
these absorptions, located at 513.4, 535.3, 556.3, 583.0, 604.9, 
and 639.2 nm, are assigned to the 2 A"(0, v", 0) «- X 2 A'(0,0,0) 
H CO' (y"=8-13) transitions ba sed on band positions reported 
by Ivan IJzendoorn et al.1 dl983l) . The clear HCO progression 
indicates a photolytic source of free H atoms in the ice. In ad- 
dition, it confirms the ability of this setup to record reactive 
intermediates in the ice. 

A possible explanation of the source of H atoms is re- 
lated to the experimental conditions. The experiments reported 
here are performed under high vacuum (10 7 mbar) condi- 
tions. Therefore, background H2O vapor has ample time to 
condense onto the sample window while cooling down and 
while growing the ice sample. Water is well known to pho- 
todissociate upon VUV irradiation (e.g. IQberg et al. 2009; 
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Fig. 3. Vibrational progression of HCO' generated in a Py:CO 
ice at 10 K after 600 seconds of VUV irradiation plotted to- 
gether with a Gaussian fit (red) to the absorption spectrum 
(black). The individual Gaussians are shown in blue. The or- 
ange line indicates the red wing of the underlying broad ab- 
sorption feature (section [3~6l l. (This figure is available in color 
in electronic form.) 



Andersson & van Dishoeckll2008h : 



VUV 

H 2 Q > OH +H. 



(3) 



An experiment on VUV irradiation of a "pure" CO ice indi- 
cated that HCO is also efficiently produced in the absence of 
pyrene. Therefore, it is likely that water contamination is re- 
sponsible for the production of HCO' via: 



H + CO -> HCO . 



(4) 



Another possible formation route could be via VUV induced 
hydrogen abstraction from pyrene. This pyrene photodissoci- 
ation reaction, however, is unlikely to occur, since PAHs are 
generally highly photostable molecules. 

3.4. The 400 nm band carrier 

Another vibrational progression appears at 400.5, 392.5, and 

378.4 nm in the CO ice experiments. As shown in Fig. |U the 

400.5 nm band dominates this progression. In contrast, a sin- 
gle band appears at 399.4 nm in the Py:H20 ice upon VUV 
irradiation of the samples. The relative intensity of these bands 
varies with respect to the Py' + bands. The 400 nm bands are 
more pronounced than the cation bands in the H2O ice only at 
high temperatures, whereas they are more pronounced in the 
low temperature CO ice. 

Two additional measurements were performed to identify 
the carrier responsible for these transitions. A kinetic experi- 
ment was performed on non-VUV irradiated Py:CO ice. This 
ice showed no sign of pyrene ionisation by the Xe-lamp which 
is used as spectroscopic light source, nor of the production of 
HCO' and of formation of the 400 nm band. Subsequently, the 
ice was irradiated with the VUV source for 10 minutes. The 





-5 -10 -15 

Integrated absorbance Py 



Fig. 4. Integrated absorbance of the 400 nm PyH band growth 
plotted against the loss of the 334 nm Py band after VUV irra- 
diation is stopped in a at 10 K ice. The straight line directly re- 
flects the relative oscillator strength of both bands as described 
in Section l3~2l The two insets show the PyH and Py vibrational 
progressions in a CO ice 90 minutes after photolysis is stopped. 



steady growth of the 400 nm band with VUV photolysis in- 
dicates that the species responsible for the 400 nm band is a 
product of VUV processing of the ice. Moreover, when the 
VUV irradiation is stopped, the 400 nm band carrier contin- 
ues to grow at the expense of the remaining neutral pyrene. 
This indicates that the chemical reaction leading to the forma- 
tion of this species is not directly photon dependent, but rather 
depends on the diffusion of a photoproduct. A similar experi- 
ment on a Py:H20 ice at 25 K indicates that the same process 
also takes place in H2O ice. The detection of HCO radicals in 
the ice and the inherent presence of free photolytic H atoms, 
hints that the growth of the vibrational progression starting at 
~400 nm could be the result of the reaction of pyrene with dif- 
fusing H atoms: 

Py + H -> PyH. (5) 

This assignment to the 1 -hydro- 1-pyr enyl radical (PyH '^ is sup 

ll 



ported by other experimental studies ( Okadaetal Jl97fill980h . 
where progressions at similar band positions are observed upon 
(laser) flash photolysis. 

Contrary to the Py:H20 experiments where pyrene is also 
efficiently ionized, the experiment on PyH formation in CO 
shows no sign of other reaction products. The integrated ab- 
sorbance of the growing PyH' transition is plotted versus the 
integrated absorbance of the diminishing neutral in Fig. @] 
Growth is tracked over a duration of more than 1 .5 hours. Since 
there is a one-to-one conversion of Py to PyH' in the Py:CO 
ice (Eq.|5]l, as described in Section [3721 we derive an oscillator 
strength of 0.089 for this species by fitting a line through the 
correlating absorbances in Fig. [4] 

3.5. The 405 nm band carrier 

Besides the Py + and the PyH bands, another distinct absorp- 
tion is found in the spectra of VUV irradiated ices. This feature 
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is located at 405.0 nm in the Py:H20 and at 406.2 nm in the 
Py:CO experiment. In our previously published paper on low 
temperature Py:H20 ice we tentativel y assigned this absorp - 
tion to a negative ion, Py~ or PyO (IBouwman et al.l 120091) . 
The experiments on Py:CO ices presented here enable us to 
rule this assignment out because of the nearly absent Py' + tran- 
sitions. Firstly, Py" is ruled out because a muc h stronger sec- 
ond P y~ absorption band, expected at 490 nm dMontejano et al.l 
19951) . is absent in our Py:CO experiment. Secondly, PyO is 
also rule d out, because it should exhib it absorptions down to 
350 nm (Milo savlievic & Thomasll2002 ). bands which are also 
absent in the Py:CO experiment. Additionally, in our previous 
paper we assumed that PyCT was a product of PyOH. The for- 
mation of PyCT is also unlikely in the absence of PyOH ab- 
sorptions in these experiments, as discussed below. 

The absorption at 405 nm does not correlate with the cation, 
nor with the PyH' band. The band only appears during pho- 
tolysis and hence is characterized as a VUV photon related 
product. From the literature it is known that a pyrene triplet- 
triplet ( 3 A ? <— 3 -Bj M ) transition is expected at this wavelength 
upon laser excitation of pyrene in solution w hich populates the 
lowest member of the triple t manifold (e.g. iHsiao & Webber! 
1992 ; ILangelaar etaL I ll970h . In order for the 405 nm band to 
arise from this triplet-triplet transition, the lowest level must be 
populated and remain so with a long enough lifetime to allow 
absorptions to the 3 A~ level. In the ice experiments reported 
here, there are a number of possible routes for pumping the 3 Py 
state. The most obvious route is via a photo excitation followed 
by intersystem crossing: 



!p y 



VUV 



'Py* 



3 Py. 



(6) 



Triplet formation is fou nd to decrease with decreasing temper- 
ature in an ethanol ice (IStevens et al.l l 1967). This translates to 
our experiment in a nearly absent 405 nm band in the low tem- 
perature Py:H20 experiment, because of the high Py + forma- 
tion efficiency. In the high temperature H2O ice experiments, 
on the other hand, the 405 nm absorption is much stronger be- 
cause pyrene is available. 

In the CO ice, on the other hand, where Py + production is 
low, formation of the 405 nm band carrier appears to be very 
efficient at low temperatures. The production of the 405 nm 
band carrier needs VUV photons to be initiated. Again, pump- 
ing of the 3 Py state can occur via Eq. [6] Moreover, CO has a 
dipole allowed electronic transition in the VUV. Hence, specu- 
lating, pumping of the 3 Py state by collisional de-excitation of 
CO molecules exited by the VUV radiation provides a reaction 
path: 



followed by: 



VUV 

CO > CO* 



CO* + 1 Py -> 3 Py + CO. 



(7) 



(8) 



In summary, while we cannot identify the carrier of the 

+ 

2// 



3.6. Broad absorption feature 

Finally, we discuss a broad feature extending from about 350 
to 570 nm that underlies the narrower bands reported in the 
previous sections that grows upon photolysis in all cases. This 
band probably arises from overlapping bands due to a num- 
ber of Py/H20 or Py/CO photoproducts. Part of this Py-residue 
feature remains even after warming up the sample window to 
room temperature, whereas all other features disappear at the 
water desorption temperature. 

As discussed above, the very broad feature must be due to a 
variety of similar but distinct photoproducts, all containing the 
pyrene chromophore. Mass spectral analysis of the species pro- 
duced by the VUV photolysis of a few other PAHs in water ice 
show that the parent PAH is not destroyed but t hat OH, O, and 



H are added to some of the edge carbon atoms dBernstein et al 



405 nm band, the 3 A <— ^Bi transition seems a plausible ex 



planation. 



19991) . Given the multiplicity of the side sites on pyrene that 
can undergo substitution, it is likely that the photoproducts pro- 
duced in the experiments reported here are multiply substituted, 
rather than singly substituted. Thus, it is plausible that a mix- 
ture of related but distinct Py-X„ species, where X may be H, 
OH, or O, produce the broad band. 

In our previous work we reported the production of a clear 
and reproducib le PyOH band at 344. 9 nm in a low temper- 
ature H2O ice dBouwman et alj 2009). The results presented 
here do not show this absorption feature. However, in some 
instances the absorption was found upon irradiation or warm- 
up of the ice. The irregular appearance of the PyOH absorp- 
tion feature in these experiments indicates that the formation of 
this species is highly sensitive to the sample's physical param- 
eters, i.e. structure of the ice, temperature and concentration. 
One possible explanation is that in the previously reported ex- 
periments, the Py concentration was not controlled and those 
experiments sampled a very different ice concentration and, by 
implication, physical ice structure. While we do not have a so- 
lution for this discrepancy, we would like to stress that both 
measurement series have been fully reproducible over many in- 
dependent experiments for periods of months. An experimental 
program to investigate the role of concentration on PAH:H20 
ice photochemistry is underway. 

4. Py:H 2 ice photochemistry at different 
temperatures 

Figure [5] shows the spectral evolution of two different Py:H20 
samples at different temperatures. The top frame presents the 
280 to 540 nm spectra of the 100 K Py:H 2 ice after 0, 20, 40, 
80 and 160 s of in situ photolysis and the bottom frame the cor- 
responding spectra for the 25 K ice. These spectra are snapshots 
from the more than one thousand spectra collected during 4 hr 
of photolysis. They illustrate the rapid changes that take place 
during the early stages in the photochemistry of these ices and 
the major differences in reaction products at different tempera- 
tures. 

To probe the VUV driven photophysics and reaction dy- 
namics for a set of selected temperatures, the production and 
depletion of species was tracked as a function of irradiation 
time. To this end, the Py 334 nm, Py' + 445 nm, and the PyH 
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Fig. 5. The VUV induced spectroscopic changes in Py:H20 ice 
for two different temperatures as a function of photolysis time. 
Comparing the spectra from the 25 K ice (bottom) with those 
of the 100 K ice (top) shows the critical role that temperature 
plays in determining photochemical pathways in a PAH con- 
taining ice. In the 25 K ice, cation formation is favored over 
production of the pyrene residue and the 400 and 405 nm band 
carriers. The opposite holds for the 100 K ice. 

400 nm bands were integrated for every spectrum. The spec- 
tra in Figs. [1] and [5] show that it is rather straightforward to 
determine the boundaries needed to integrate these bands. We 
estimate that the uncertainty in most of these band areas is of 
the order of 10%. 

The integrated absorbances of the neutral Py, strongest Py' + 
and PyH bands in H2O ice at temperatures of 25, 50, 75, 100 
and 125 K are plotted versus photolysis time (VUV fluence) 
in Fig. [6] The spectra in Fig. [5] and photochemical behavior in 
Fig.|6]show that, upon photolysis, neutral pyrene loss is imme- 
diate and rapid. The initial growth of Py + mirrors the rapid, ini- 
tial loss of Py. However, while Py steadily decreases, and sev- 
eral other Py photoproduct bands increase during some 4 hours 
of photolysis, the production of Py + reaches a maximum and 
then slowly diminishes. From Fig. [6] one can clearly see that 
ionisation of pyrene is most efficient in the low temperature ice. 
Formation of PyH', on the other hand, is much more efficient at 
higher temperatures. 

For comparison, the integrated absorbances for the irradi- 
ated Py:CO ice are plotted as a function of time in the right 
bottom frame of Fig. [6] Note that the PyH' band is multiplied 
by a factor of 10 in the Py:CO experiment, compared to a fac- 



tor of 20 for the Py:H20 experiment. The PyH' band is clearly 
more prominent in the CO ice experiment compared to the H2O 
ice experiments. The Py + signal on the other hand is negligible. 
This indicates that the H2O ice plays a role in ion formation and 
stabilization. 

To place this behavior on a quantitative footing, the inte- 
grated areas for the Py and Py ' + bands are converted to number 
densities using Eq. Q] Here, an oscillator strength of 0.33 is 
used for the 334 nm Py bands. The values used for the oscil- 
lator strengths of the Py + and PyH bands are 0.33 and 0.089, 
respectively, as determined in Sections |3 .21 and [3~4l Perusal of 
Fig.|6]shows that Py behaves similarly in all of the ices consid- 
ered here. Regardless of temperature, its signal drops quickly 
with the onset of irradiation and continues to diminish with 
ongoing photolysis. Likewise, Py + grows rapidly with initial 
photolysis but it peaks after a relatively short time interval cor- 
responding to a fluence of roughly 8 x 10 16 photons and then 
drops continuously. While the Py' + growth and loss curves re- 
semble one another, cation production efficiency is strongest 
in the 25 K ice. This efficiency remains of the same order at 
even lower temperatures (not shown here). The photolysis time 
required for the cation to peak shortens with temperature in- 
crease. The PyH' band contribution is minor with respect to the 
Py' + band for ices below 50 K. This reverses between 50 and 
75 K, suggesting that there is a change in the dominant Py:H20 
ice photochemical channel in this temperature range. 

A kinetic analysis of the plots in Fig.|6]is carried out using 
the reaction scheme indicated in Fig. [7] Here, k\\ is the pho- 
toionisation rate of Py to Py ' + , kn the electron-ion recombina- 
tion rate of Py + , £21 the production rate of the PyH' feature, and 
£22 the reverse reaction rate of PyH' to Py. The rates designated 
k\, k%, and £3 are the production rates for the different products 
that comprise the Py-residue band. The oscillator strengths for 
the Py + and PyH' bands are also fitted, but are restricted to 
stay within +10% of the experimentally determined values of 
0.33 and 0.089. All reactions are assumed to be first order in 
the reactant. The relative abundances of "free or solvated elec- 
trons" and O, H and OH radicals in the ice are not considered. 

The fits to the growth and decay curves are included in 
Fig. [6] and the temperature dependence of the derived rate con- 
stants is presented in Fig. [8] The agreement between the fit and 
the experimental data in terms of curve shape and absolute in- 
tensity is good. The fitted oscillator strengths of the Py' + and 
PyH' bands amount to 0.31 and 0.082, respectively, and hence 
do not deviate much from the experimentally determined val- 
ues. 

The graph in Fig. [8] indicates that the Py photoionisa- 
tion rate (k\\) drops rapidly between 25 and 50 K. The elec- 
tron recombination rate (k\i) only decreases slightly, if at all, 
within the error over the entire temperature range. As men- 
tioned above, the production of the PyH becomes more im- 
portant at higher temperatures. Its formation rate (&21) is small 
in all ices up to 50 K (< 4.4 x 10~ 5 ), but jumps to > lxl0~ 4 
in the ices with temperatures of 75 K and higher. Also the back 
channel from PyH' to Py, £22, shows a temperature dependence. 
It increases nearly linearly in going from cold to warm ices. 
The formation rate of a photoproduct directly from Py {k\) also 
seems to jump at 50 K. The formation rate of products from the 
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Fig. 6. The integrated absorbance of the Py 334 nm, Py + 445 nm, and PyH 400 nm bands as a function of VUV irradiation 
in Py:H20 ices at 25, 50, 75, 100, and 125 K and a Py:CO ice at 10 K plotted together with the fits (grey lines) described 
in Section |4] Integrated absorbance values are scaled and normalized to the initial value for the Py signal. The PyH' feature is 
multiplied by a factor of 20 for the Py:H20 experiments and by a factor of 10 for the Py:CO experiment. The approximate, overall 
growth of the total Py photoproduct band (P1+P2+P3) is also shown (dotted line). (This figure is available in color in electronic 
form.) 



Py + species, on the other hand, seems to lower with increas- 
ing temperature. Finally, the rate of product formation from the 
PyH channel is low throughout the entire temperature range. 
The jump in rate of the formation of Pj and PyH with temper- 
ature probably reflects the diffusion barrier of radical species 
(H and OH) in the ice. 

Since published studies of the processes induced by the 
photolysis of other PAH:H20 ices are limited, there is not 
much information available with which to compare these re- 
sults. While, to the best of our knowledge, there are no reports 
on the photochemistry that takes place as a function of ice tem- 
perature or long term fluence, the VUV photochemistry of the 
PAHs naphthalene, 4-methylpyrene (4MP), and quatterrylene 
in water ice at 10 K has been studi ed (Gudipati & Allama ndola 
20031: lGudipatil2004l:lGudipafi & A llamandola 2006aH). Their 



results are in good agreement with the low temperature (25 K) 
case reported here. Namely, the parent PAH is easily and ef- 
ficiently ionized, with quantitative conversion of the neutral 
species to the cation form. The focus in the earlier studies was 
on cation production and stabilization and not on long duration 
photolysis ex periments. In their study of 4 MP:H?Q (1:>500) 
ice at 15 K, iGudipati & Allaman dola (2003) utilized a reac- 
tion scheme similar to that on the right half of that presented 
in Fig. [7] Table [2] compares the reaction rates they determined 
with those of the 25 K ice reported here. Except for the pro- 
duction of P2, that deviates by one order of magnitude, there is 
very good agreement between the rate constants for each step 
in the two experiments. 

The growth and decay curves in Fig. [3 taken together with 
the temperature dependence of the reaction rates in Fig. [8] show 
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Fig. 7. Reaction scheme used to fit the experimental data. 
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Fig. 8. Parameters (k„ m ) as a function of temperature resulting 
from fitting the reaction scheme (Fig. |7]i to the kinetic experi- 
ments (Fig. [5). All rates are indicated in s . 



that the VUV driven PAH photochemistry depends strongly on 
ice temperature. The influence of the ice morphology on this 
chemistry was also investigated, to understand the origin of the 
temperature dependence. An experiment on an ice deposited 
at 25 K, annealed to 125 K, and subsequently cooled down to 
25 K before photolysis, showed that the ionisation rate and effi- 
ciency are similar to that of an unannealed ice. Apparently, it is 
not the morphology but the temperature of the ice that primarily 
determines which process dominates. We discriminate between 
two temperature regimes. One governed by ion-mediated pro- 
cesses that dominate at 25 K and slightly higher temperatures 
and a second, presumably radical driven regime, that becomes 
increasingly more important at higher temperatures. 

5. Astrochemical Implications 

As shown in the previous sections, ionisation and chemistry 
of a rather small PAH, pyrene, trapped in H2O ice turns out 
to be very efficient in a laboratory setting. Here, we extend 
these findings to interstellar conditions, with the aim to include 
rates in astrochemical models. For this, it is crucial to disen- 
tangle pure photochemical processes from diffusion, since the 



Table 2. The reaction rates for the VUV photoly- 
sis of Py:H 2 (-1:5,000) ice at 25 K compared to 
those for 4-mefhylypyrene:H 7Q (1:<500) ice at 15 K 
(Gu dipati & Allamandolal2003l) . Photon rates are indicated for 
reaction channels which are dominated by photon processes. 



Rate 




This work 


Photon rate 


Gudipati 2003 






(s- 1 ) 


(crrrphotorT 1 ) 




VUV 

Py >Py + 


(*n) 


(1.2±0.1)xl0 -3 


1.2x10"" 


1.3xl0" 3 


Py + +e" -^Py 




(9±2)xl(r 4 


9xl0" 18 


8xl0- 4 


Py -> Pi 


(Jfcl) 


(5±l)xl0" 5 




4x1 0~ 5 


Py + -> P 2 


ft) 


(5±l)xl(T 4 




5xl0" 5 



latter will be highly dependent on the number density of radi- 
cals and electrons in the ice. As mentioned in the previous sec- 
tion, the photoionisation of Py is probably a single photon pro- 
cess, whereas protonation of Py and the electron recombination 
of Py' + are results of both VUV photolysis and diffusion. The 
mechanism for PyH deprotonation is not clear, since it can ei- 
ther proceed through VUV processing or through hydrogen ab- 
straction by diffusing species. Diffusion of radicals through the 
ice is a thermally activated process and will therefore increase 
with temperature. Recombination, however, is largely temper- 
ature independent in our experiments, indicating that the rate 
of Py + recombination is not dominated by diffusion of elec- 
trons in the ice. If Py' + loss is via electron recombination and 
not Py + reaction with H2O or one of its photoproducts, the 
electron most likely originates from the initial photoionisation 
event where electrons stay in the vicinity of the recombining 
Py' + species. Hence, this local process can be, although indi- 
rectly, regarded as a single photon process. 

The rates of protonation of Py and deprotonation of PyH 
show a temperature dependence and the importance of diffu- 
sion can therefore not be excluded. This makes it harder to di- 
rectly translate the rates (s _1 ) to photon rates (cm 2 photon _1 ). 
However, we can determine astrochemical photon rates for 
both ionisation and recombination of pyrene in interstellar H2O 
ice (see Tabled. 

Now, to translate this to the astrochemical situation and 
compare this to other processes, let us assume that PAHs gen- 
erally have ionisation rates similar to that of pyrene. How do 
ionisation and chemistry compare with other processes such as 
photodesorption of the icy grain mantle, in which the PAHs are 
embedded? To exemplify this, the rate of ionisation of a PAH 
in water ice at 25 K (in photon ') is calculated anywhere in 
a dense cloud where A v = 3 and com pared to the VUV pho- 
todesorption rate of H2O as derived bv lOberg et al. I d2009h . It 
is well established that the onset of ice formation occurs in 
clouds with an edge-to -edge (through the cloud) magnitude 
of A v = 3 (e.g. IWhittet et alJl200lh . Thus inside our hypo- 
thetical dense cloud at A v - 3 (from cloud edge to within 
the cloud) ices are present. 

The experimentally determined PAH ionisation rate in H 2 
at 25 K, normalized to the total amount of deposited PAH is 
given by: 



kn = 



A [PAH + ] 
[PAH]p 

dr 



= 10 3 s" 1 . 



(9) 
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Consider a typical interstellar grain, covered by a 100 mono- 
layer (ML) thick ice. The number of sites on a grain is 
10 15 cirT 2 . If we assume that one in every 10 4 particles on the 
grain is a PAH, the total number of PAH molecules on the grain 
is [PAH] = 100 • 10 15 • 10" 4 = 10 13 cirT 2 . Furthermore, the 
VUV photon flux in our laboratory, <D, is 10 14 photons cm"^" 1 . 
The production rate of PAH cations on an interstellar grain is 
now given by [PAH]o • = 10" 4 photon -1 . This ionisation 

rate is an order of magnitud e lower than the rate of photodes- 
orption (~ 10 3 photon" 1 ) dQberg et al.ll2009h . 

However, in our dense cloud the number of photons avail- 
able for PAH photoionisation is larger than the number of 
photons available for photodesorption of H2O ice. This is be- 
cause H2O photodesorption occurs primarily by VUV photons, 
whereas PAH ionisation can occur with much lower energy 
photons. To quantify the radiation field in a dense cloud at 
Ay = 3 as a function of wavelength (A ), we take t he ave rage 
UV interstellar radiation field (I v ) from [Sternberg (1988) and 



rewrite the expression to 1^ with units photons cm s nm 
1.068 x 10" 4 c 1.719 xl0" 2 c 6.853 x 10"^ 



A? 



A* 



A 5 



(10) 



where c is the speed of light in nm s 1 . The attenuation of the 
radiation field by dust as a function of wavelength is given by: 



-U.6A V A A 
D A = exp( ), 



(ID 



"V ^lOOOA 

from lDraine & Bertoldil 1 1996), whe re we assume R v = 3.1 
and Av/A mok = 0.21 dWhittetll2003l) . This results in: 



D A =exp(-0.8— -A v ), 
Av 



(12) 



where the table of A A /Av values is taken from Mathis (1990). 
The photon flux per second per wavelength interval is given by: 



(13) 



Water ice absorbs p hotons with wa velengt hs ranging from 



130 to 150 nm (Kobavashi 1983; Andersson & van Dishoeck 



2008). The ionisation energy of PAHs on the other hand, is low 
ered by about 2 eV when i n H 2 ice dGudipati & Allamandola 
20041 IWoon & Parkll2004 . For the wavelength range available 
for ionisation of PAHs, assuming that H2 O blocks all photon s 
below 150 nm, we take 150 to 250 nm dLi & Drainell200lb . 
Integrating the photon flux in a cloud of Ay = 3 over both 
wavelength intervals gives a number of photons available for 
PAH ionisation which is 6 times larger than the number of pho- 
tons available for photodesorption of H2O . Additionally, at 
Av = 3, the cosmic ray induced UV field is negligible com- 
pared to the interstellar UV field flShen et al.ll2004l) . Therefore, 
the occurrence of photoionisation is of similar order as pho- 
todesorption of the main component in the grain mantle in a 
dense cloud. The ionisation rates from Table [2] can be directly 
included in astrochemical models in the form: 



d[PAH + ] 
dr 



= fcii¥[PAH], 



(14) 



where [PAH + ] is the concentration of the PAH (pyrene) cation 
in the ice, k\ \ is the photon rate in cm 2 photon" 1 , *P is the photon 



flux in photon s^cm 2 , and [PAH] is the PAH (pyrene) concen- 
tration in the ice. 

In the above calculation we assume that all PAHs exhibit 
the ionisation behavior of the pyrene chromophore. Of course, 
more PAHs need to be investigated experimentally in order to 
draw conclusions on their general photochemical behavior in 
interstellar ices. However, if all PAHs have ionisation rates sim- 
ilar to pyrene, photoionisation and subsequent chemical reac- 
tions of PAHs trapped in ices are important processes in dense 
clouds. When frozen out in ices, PAHs have an important im- 
pact on the radical and electron budget in solid state chemistry. 
Hence, the processes described here may be more important 
than previously assumed in modeling complex interstellar grain 
chemistry. 

6. Conclusions 

A recently constructed setup has been used to track, on a sub- 
second timescale, the photochemistry of a PAH in H2O and CO 
ices as a function of temperature. The setup used here clearly 
has advantages compared to the relatively slow infrared pho- 
tochemical ice studies. The conclusions from this work on the 
PAH pyrene trapped in an interstellar ice analogue are summa- 
rized below: 

1. A set of photochemical reaction products has been as- 
signed, both in irradiated Py:H20 and Py:CO ice experi- 
ments. The reaction products result from direct photoioni- 
sation of pyrene, or from a reaction of the parent, pyrene, 
with free H atoms produced in the matrix. Additionally, an 
absorption band is tentatively assigned to a triplet-triplet 
transition of pyrene. A vibrational progression assigned to 
HCO' is found in spectra of the VUV irradiated Py:CO ice. 

2. Pyrene is easily and efficiently ionised when trapped in 
H2O ice. Photoionisation is a non-diffusion related reaction 
and hence a photon-rate of 1.2x 10" 17 cm 2 photon" 1 , which 
can serve as input for astrochemical models, is derived. 

3. When trapped in CO ice, pyrene ionisation is inefficient 
compared to water ice. 

4. Electron-ion recombination is independent of ice tem- 
perature and is characterized as a non diffusion domi- 
nated reaction. For this process, a photon rate of 9 x 
10" 18 cm 2 photon" 1 is derived which can be directly used 
in astrochemical models. 

5. There are two distinct reaction paths in the photochemistry 
of pyrene trapped in H2O ice. At low temperatures (< 50 K) 
the chemistry is dominated by ion-molecule interactions 
and processes. At temperatures above 50 K reactions are 
dominated by diffusing radical species. 

6. A simple model indicates that, in dense clouds where 
Av = 3, the rate of pyrene ionisation is comparable to the 
rate of photodesorption in water rich ices. Hence, chemical 
reactions involving pyrene, and its cation can be important 
in modeling grain chemistry in such environments. 
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